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摘 要 :在 渤海 和 北 黄海 海域 , 冰 和 载荷 是 影响 风电 基础 结构 振动 响应 及 疲劳 寿命 的 关键 因素 。 对 于 
单 桩 风电 柔性 结构 而 言 , 在 水 线 处 安装 抗 冰 锥 可 有 效 降低 冰 载 荷 , 保 障 海 上 风电 安全 运行 。 为 分 析 
海 冰 对 锥 体 风 电 结构 的 冰 载 荷 , 基 于 扩展 多 面体 离散 元 方法 模拟 了 平整 冰 和 碎 冰 对 锥 体 风电 结构 
的 作用 过 程 ,分 析 了 不 同 锥 径 和 海 冰 密集 度 下 的 锥 体 风 电 结 构 冰 载荷 及 冰 激 结 构 振 动 响应 。 计 算 
结果 表明 , 锥 体 结 构 冰 载荷 具有 显著 的 随机 性 和 周期 性 ,平整 冰 对 锥 体 结 构 冰 载荷 要 明显 高 于 碎 
冰 ; 平 整 冰 对 上 锥 体 的 冰 载 荷 大 于 下 锥 体 冰 载荷 ,而 碎 冰 对 上 锥 的 冰 载 荷 却 小 于 下 锥 体 冰 载荷 ;所 
5 提 方法 可 有 效 模拟 海 冰 对 锥 体 风 电 结 构 的 作用 力 。 
关键 词 :离散 元 方法 ; 冰 载 荷 ; 锥 体 风 电 结构 ;平整 冰 ; 碎 冰 区 
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Discrete element analysis of ice loads on conical 
monopile structure of wind turbine in level ice 
and broken ice field 
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Abstract:The ice load is a key factor affecting the vibration response and fatigue life of wind turbine 
structure in the Bohai Sea and the North Huanghai Sea of China. For the flexible monopile of wind tur- 
bine ,the anti-ice cones installed on the structure at the water level can effectively reduce the ice load, 
which can ensure safe operation. To analyze the ice load of conical pile of wind turbine ,the dilated poly- 
hedron based discrete element method ( DPDEM) is adopted to model the level ice and broken ice floes. 
The ice loads and ice-induced vibration of conical structure are simulated with various conical diameters 
and ice concentrations of broken ice field. The simulated results show that the ice loads of level ice and 
broken ice perform the characteristics of randomness and periodicity. The ice load of level ice is obvious 


higher than that of ice floes. The ice load of level ice on the upward cone is larger than that on the down- 
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ward cone, while the ice load of broken ice on the upward cone is smaller than that on the downward cone. 


The study shows that the DPDEM can be applied to simulate the ice load on conical structure of wind tur- 


bine effectively. It can be applied to the ice-resistant structure design and safe operation of offshore coni- 


cal structures in ice-covered regions. 


Key words discrete element method ; ice load ; conical structure of wind turbine foundation ;level ice ; bro- 


ken ice field 


随 着 可 再 生 能 源 的 迅速 发 展 ,风力 发 电 已 成 为 
其 不 可 或 缺 的 重要 组 成 部 分 。 相 较 于 陆地 风力 发 
电 , 海 上 风能 因 储备 资源 更 为 丰富 平均 空气 密度 高 
且 不 占用 土地 资源 等 优点 而 具有 广阔 发 展 前 景 "1。 
在 渤海 及 北 黄海 , 冬 春季 节 受 西伯 利 亚 强 冷 空气 影 
响 , 使 得 风能 密度 远 高 于 其 他 海域 。 因 此 , 寒 区 风力 
发 电 成 为 我 国 未 来 发 展 海上 风电 的 重点 。 在 高 纬 
度 守 冷 海域 , 海 冰 不 可 避免 地 与 风电 基础 结构 发 生 
相 王 作用 ,产生 强烈 的 冰 载 荷 和 结构 动力 响应 = 。 
单 往 风 机 结构 作为 一 种 柔性 细 长 结构 , 相 较 于 海上 
石 消 平 台 等 传统 海洋 工程 结构 ,其 在 海 冰 作用 下 的 
并 稚 响 应 更 为 复杂 。 在 随机 振动 下 ,风机 结构 的 
冰 小 荷 及 动力 响应 对 冰 厚 和 海 冰 强度 较为 敏感 , 严 
恒 影 响 风 机 结构 的 疲劳 寿命 与 运行 安全 9 。2007 
^s Huang 等 "通过 室内 冰 池 实验 开展 了 不 同 冰 速 
下 的 单 桩 结构 冰 激 振动 研究 。 
寒冷 海域 的 风电 基础 结构 需要 具备 良好 的 抗 冰 
性 能 。 由 于 海 冰 的 压缩 强度 远大 于 其 弯曲 强度 , 锥 
体 结构 能 有 效 将 海 冰 挤 压 破坏 转化 为 弯曲 破坏 ,从 
而 恬 低 结构 冰 和 载荷 和 冰 激 结构 振动 中 。 风 机 结构 的 
冰 材 荷 不 仅 与 冰 厚 、 冰 强度 、 冰 速 有 关 , 同样 与 海 冰 
分 所 类 型 锥 径 等 参数 密切 相关 。 龙 雪 等 通过 球 
体 离散 元 模拟 发 现下 锥 结构 可 有 效 降低 海洋 平台 的 
冰 载 荷 。2013 年 , Tian 等 ' " 通过 室内 模型 试验 发 
现 , 当 锥 径 与 平整 冰 厚 之 比 小 于 30 时 ,其 冰 载 荷 随 
锥 径 的 增 大 而 增 大 。 而 有 关 平 整 冰 与 碎 冰 区 对 风电 
锥 体 结构 冰 载 荷 相互 作用 的 比较 研究 较 少 。 

离散 单元 法 5 广泛 应 用 于 岩 土工 程 、 化 工 工 
程 和 海洋 工程 领域 ,特别 是 在 寒 区 海洋 工程 方面 。 
离散 元 方法 因 其 能 够 有 效 描述 海 冰 破碎 运动 过 程 中 
的 动 冰 载 荷 而 被 广泛 应 用 ,并 发 展 了 球体 ”、 赂 
RUU 多 面体 "和 扩展 多 面体 "等 不 同 单元 类 型 。 
基于 球体 单元 的 黏 结 -破碎 模型 具有 模型 简单 .计算 
效率 高 等 优点 ,可 模拟 海 冰 与 海洋 工程 结构 的 相互 
作用 与 破碎 过 程 ”” ;基于 多 面体 模型 离散 元 方法 
可 有 效 描述 平整 冰 对 海洋 工程 结构 相互 作用 动力 过 


程 中 的 弯曲 破坏 、 堆 积 重 琶 现 象 ; 基于 扩展 多 面 
体 模型 的 离散 元 方法 能 够 很 好 地 描述 海 冰 与 海洋 工 
程 结 构 相 互 作用 的 动力 过 程 ,并 且 可 提高 多 面体 单 
元 与 海洋 结构 作用 的 计算 效率 ,在 海 冰 离 散 元 领域 
具有 明显 优势 。 

本 文采 用 扩展 多 面体 离散 元 方法 分 别 构造 平整 
冰 及 碎 冰 区 ,模拟 海 冰 与 锥 体 风电 基础 结构 的 相互 
作用 过 程 ,并 将 计算 冰 载 荷 结果 与 EC 规范 和 ISO 标 
准 进行 对 比 ,以 验证 其 可 靠 性 ,进而 研究 平整 冰 及 碎 冰 
对 锥 体 风电 基础 结构 的 振动 响应 及 冰 载 集 变 化 规律 。 


1 扩展 多 面体 离散 元 方法 


采用 扩展 多 面体 离散 元 方法 可 以 有 效 模拟 海 冰 
的 几何 形态 ,并 通过 单元 间 的 秋 结 -破碎 模型 描述 海 
冰 与 单 桩 锥 体 风电 基础 结构 的 破坏 过 程 ,由 此 确定 
不 同 冰 况 和 结构 参数 影响 下 的 冰 载 集 。 


1.1 海 冰 扩展 多 面体 离散 元 方法 


若 任意 给 出 空间 中 的 两 个 几何 体 点 集 4 和 B, 那 
么 空间 集合 4 和 B 的 Minkowski Sum 定义 为 
AQB = ix *yIx CA,y CB] (1) 
式 中 ,x 和 y 分 别 为 集合 4 和 B 中 的 几何 坐标 点 。 
空间 体 4 M B 分别 定义 为 基本 多 面体 与 扩展 球 
体 , 基 本 多 面体 与 球体 单元 经 过 Minkowski Sum Æ 
加 ,形成 一 个 具有 光滑 度 的 扩展 多 面体 ,如 图 1 所 
示 , 扩 展 球体 半径 称 为 扩展 半径 。 采 用 Voronoi 切割 
算法 将 一 定 区 域 的 二 维 平面 随机 划分 为 若干 任意 多 
面 形 , 其 形状 与 随机 点 分 布 位 置 密 切 相 关 "”。 将 任 
意 多 面 形 在 垂直 方向 上 拉 伸 并 基于 Minkowski Sum 
方法 构造 具有 一 定 厚度 的 扩展 多 面体 海 冰 单 元 "| 。 
由 此 生成 的 平整 冰 初 始 场 如 图 2(a) 所 示 。 对 于 碎 
冰 区 初始 场 ,可 通过 改变 海 冰 禾 盖 面 积 与 海域 总 面 
积 之 比 即 海 冰 密 集 度 ,生成 具有 随机 形态 和 面积 分 
布 的 碎 冰 区 ,如 图 2(b) 所 示 。 
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图 1 基于 Minkowski Sum 的 扩展 多 面体 
Fig.1 Dilated polyhedral element based on Minkowski Sum 


(a) 平 整 冰 初始 分 布 


(b) 碎 冰 区 初始 分 布 
图 2 基于 扩展 多 面体 离散 元 的 海 冰 初始 分 布 


Fig.2 Initial distributions of level ice and broken ice field 


constructed with the dilated polyhedral discrete elements 
二 在 扩展 多 面体 离散 元 的 接触 求解 过 程 中 ,可 采 
用 基于 近似 包 络 函 数 的 搜索 判断 算法 计算 接触 点 与 

量 ,并 采用 非 线性 的 接触 力 模型 简化 计算 扩展 
多 面体 间 的 接触 力 以 提高 计算 效率 。 法 向 接触 力 可 
dor Ari ey B TUR E 71 2 RU , 即 


& a t EA. G 
N 本 二 
AE* /R 
E k, = 3 ,Cn 76, v mak, , 
~ _ “he — 
SS 0 74 
Xon. 为 法 向 接触 刚度 ;6, 为 两 个 多 面体 单元 间 


Bp e] BH, 为 法 向 阻尼 系数 ;6, 为 法 向 相对 
HR E 和 及" 分 别 为 等 效 旨 性 模 量 和 等 效 闫 入 半 
B e, 为 量 纲 为 一 的 阻尼 系数 ;ms 为 等 效 质量 ;e 
为 单元 的 回 弹 系 数 1。 

切 向 接触 力 与 切 向 重 闯 量 和 摩擦 力 有 关 , 一 般 
最 大 切 向 力 不 超 过 摩擦 力 。 切 向 接触 力 TR 
RU 


NE 
F,-Fi «Feud Fr - (1o) | 


i A 

je AE 
式 中 : 为 切 向 弹性 力 ; F SUI BE 21 u 为 摩擦 
系数 ;6, 为 切 向 重 毒 量 ,_- 般 根据 相对 切 向 速度 并 逐 
步 大 加 的 方法 累计 计算 切 向 重 毒 量 ;5"" 为 最 大 切 
向 重 关 量 ,其 主要 由 摩擦 系数 和 法 向 重 关 量 决定 , 且 
à" =j[8,(2-v)]/(2 -2y), 这 里 为 泊 松 比 ; 


ChinaXiv 合 作 期 刊 


6, 为 切 向 相对 速度 ;c, 为 切 向 阻尼 系数 ,这 里 取 c, = 
c,/ [2(1 +v)]。 

XEFE B5 dh i i 70 , Te P] T 4H Hc ds f 7C 
间 的 交界 面 角 点 位 置 处 设 定 积分 点 。 该 积分 点 在 两 
单元 上 对 应 的 角 点 称 为 茜 结 节点 , 旦 该 交界 面 为 称 
为 黏 结 面 ,如 图 3 所 示 。 通 过 设 定单 元 的 刚度 与 阻 
尼 即 可 计算 黏 接 节点 上 的 法 向 应 变 与 切 向 应 变 。 采 
用 三 维 条 件 下 的 弹性 矩阵 , 黏 结 单元 间 的 弹性 应 力 
o Ee] 

Į% 0] E(-») . y | 
o = £= 1 -2y le 
lo Y (1 *v)(1-2v)|0 元 1 

(4) 
其 中 :5 为 弹性 模 量 ;oe 为 考虑 法 向 和 切 向 的 应 力 ， 
Hl o^ 2[o,7] ;& 为 考虑 黏 结 节 点 之 间 的 法 向 和 切 
向 应 变 , 即 e= | se, 上。 在 扩展 多 面体 单元 的 黏 结 
面 上 ,法 向 阻尼 系数 e, 与 切 向 阻尼 系数 e, RS 
作用 。 由 黏 滞 作用 产生 的 应 力 o" 和 两 个 狐 结 节点 
间 的 黏 结 力 F 分 别 写作 


eS Dei e 


4 
n 


F'-(og *g)- (6) 


式 中 : 为 应 变 率 ,上 且 e= |e,,2s| ;4 为 交界 面 面 
积 ;n 为 交界 面 上 的 医 结 节点 个 数 。 

将 求 得 的 单元 间 黏 结 力 转移 到 单元 质心 上 ,并 
考虑 流体 拖 忠 力 、 浮 力 等 其 他 作用 力 , 即 可 确定 各 离 
散 单元 质心 的 合力 。 


图 3 扩展 多 面体 单元 间 的 竹 结 面 模型 
Fig.3 The bonding surface between dilated polyhedral elements 
针对 扩展 多 面体 单元 间 的 黏 结 -破碎 模型 ,本 文 
在 法 向 拉 伸 应 力 与 切 向 剪 切 应 力 共 同 作 用 下 计算 临 
界 应 力 。 如 图 4 所 示 , 这 里 损伤 (D,) 是 指 加 载 过 程 
中 刚度 软化 的 过 程 , 即 刚度 由 大 变 为 (1-D.) 的 过 
Fe ,混合 断裂 能 是 指 在 拉 伸 和 剪 切 两 个 方向 综合 影 
响 下 的 断裂 能 G*"。 当 超出 临界 应 力 后 ,采用 损伤 度 
来 描述 材料 刚度 软化 的 过 程 ” ,并 通过 B-K 混合 由 


裂 能 模型 计算 单元 间 的 临界 变形 77, 24s nx RIT 
变形 大 于 其 临界 变形 时 , 则 该 单元 间 的 黏 结 失效 , 秋 
结 单元 发 生 分 离 。 


ô ô 


到 4 单元 间 黏 结 失效 的 应 力 和 变形 关系 


Fig.4 The stress and deformation relationship 


of bond failure between elements 
1325 单 桩 锥 体 风 电 基 础 结构 的 计算 模型 


这 里 选取 渤海 某 单 柱 锥 体 风电 基础 结构 作为 研 
完 寻 象 建立 相应 的 离散 元 模型 ,如 图 5(a) 所 示 。 本 
文 取 考虑 海 冰 对 锥 体 结构 的 冰 载 荷 , 尚 未 考虑 风 载 
荷 忆 这 里 将 风电 上 下 锥 体 交界 面 处 水 位 高 度 记 为 0 
水 富 , 取 抗 冰 锥 形 心 作为 基础 结构 振动 位 移 起 始点 。 
英信 风机 的 主要 结构 参数 为 :风机 轮 载 距 水 面 高 度 
96 记 ; 距 海面 高 度 12 m 以 上 部 分 为 薄 壁 贺 欠 位 形 结 
覆 3 风 电抗 冰 锥 位 置 处 为 薄 壁 圆 简 结构 ,直径 5.5 m. 
风纪 锥 体 结构 包括 上 锥 ,下 锥 和 抗 冰 锥 套 简 结 构 ,将 
风电 锥 体 结构 表面 划分 为 三 角形 单元 ,用 于 扩展 多 
面 傣 海 冰 单 元 与 锥 体 结构 间 的 接触 判断 。 风 电 基础 
结 税 的 主要 参数 为 :最 大 锥 径 D... = 8. 75 m; 最 小 锥 
BD ain =5.5 m; 锥 体 角度 a = 64. 9"; 锥 体高 度 刀 = 
3.95 m; 结 构 运动 刚度 =2.42 MN/m; 结构 运动 阻 
尼 c=1.67 x 105 ;结构 基 频 f=0.2826 Hz; 风 电 结 构 
总 质量 m=7.68 x 10? kg, 

为 分 析 单 桩 锥 体 风电 基础 结构 的 冰 激 振动 响 
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(a) 单 桩 锥 体 结 构 模 型 


(b) 基础 结构 运动 模型 
到 5 冰 区 单 桩 锥 体 风 电 结 构 离 散 元 模型 


Fig.5 The discrete element model of wind turbine 


structure with conical mono-pile in ice field 


2 平整 冰 区 锥 体 风 电 结 构 冰 载 谷 的 离 
散 元 模拟 及 主要 影响 因素 分 析 


考虑 潮汐 涨 落 过 程 中 ,平整 冰 对 锥 体 结构 作用 
位 置 的 差异 。 采 用 离散 元 方法 对 上 下 锥 结构 冰 载 从 
及 冰 振 响应 进行 数值 计算 和 对 比分 析 , 并 在 此 基础 
上 考虑 水 位 升降 所 对 应 的 作用 锥 径 情 况 , 讨 论 分 析 
上 下 锥 体 锥 径 影响 下 的 冰 载 荷 变化 规律 。 


2.1 平整 冰 对 锥 体 风电 结构 冰 载 荷 的 离散 元 模拟 


锥 体 风电 结构 的 冰 载 荷 与 海 冰 厚度 .速度 等 海 
冰 参 数 密切 相关 。 这 里 结合 风机 所 在 海域 的 海 冰 环 


应 ,将 基础 结构 视 为 可 做 刚体 运动 的 质量 块 ,同时 简 
化 海 冰 与 结构 相互 作用 处 的 结构 刚度 和 阻尼 ,将 结 
构 的 运动 视 为 质量 块 在 刚度 和 阻尼 影响 下 的 刚体 往 
复 运 动 ,因此 可 将 该 运动 系统 视 为 “质量 -刚度 -阻尼 
模型 ” ,如 图 5(b) 所 示 , 其 动力 学 方程 可 表示 为 

F, =mu, +cu, + ku, (7) 
式 中 :Fi 为 冰 载 位 ;m 为 结构 质量 ;ju 为 结构 位 移 ;c 
为 结构 平 动 阻 尼 ;k 为 结构 平 动 刚度 。 其 中 ,阻尼 算 
阵 取 为 Rayleigh 阻尼 形式 , 即 c = am +Bk ,参数 a 和 
B 通过 结构 的 模 态 阻尼 比 和 固有 频率 计算 , 即 a = 
0. 171 1,8 20.014 6, 


境 条 件 , 取 海 冰 厚度 有 =0.15 m, Ki v; 20.5 m/s, 
开展 海 冰 对 锥 体 结构 作用 的 离散 元 模拟 计算 。 其 主 
要 参数 包括 水 位 高 度 且 =0. 65 m、 作 用 锥 径 D = 
8.1 m, 并 参考 了 扩展 多 面体 离散 元 模拟 的 相关 文 
献 ,主要 离散 元 模拟 参数 为 : 冰 厚 h; 20.15 m; ifs 
水 流速 v=0.5 m/s; 海 水 密度 p, = 1 035 kg/m? ; Y 
冰 弹 性 模 量 E=1.0 GPa; 海 冰 摩 擦 系数 =0.2; 海 
冰 回 弹 系数 e =0. 3; 海 冰 泊 松 比 > 20. 3; 法 向 黏 结 
强度 c, = 0. 6 MPa; 海 冰 单 元 平均 尺寸 S8 20.3 m^; 
海 冰 计 算 域 1xw 246 m x 42 m; $6 82, 7] RŽ Car = 
0.6; 拖 中 力矩 系数 Ciu 20.03, 
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图 6 给 出 了 平整 冰 与 锥 体 结构 相互 作用 的 破碎 
过 程 。 从 中 可 以 发 现 ,平整 冰 呈 现 初次 断裂 、 爬 升 、 
二 次 断裂 和 清除 的 过 程 ,由 此 引起 锥 体 结构 的 交 变 
SIKER o E87 为 单 桩 风电 结构 冰 载 荷 时 程 曲线 ， 
呈现 出 很 强 的 随机 性 与 周期 脉冲 特性 。 本 文 参考 了 
文献 [27-28] 中 的 “峰值 ”提取 方法 ,该 方法 可 从 冰 
力 时 程 曲线 中 有 效 提取 “峰值 ”, 在 现场 监测 和 数值 
模拟 的 结果 分 析 中 效果 较 好 。 通 过 对 结构 冰 载 荷 时 
程 的 峰值 进行 提取 ,得 到 冰 载 集 峰 值 的 平均 值 为 
83.77 kN, 


ka 
(O (a) 现场 监测 (b) 离散 元 模拟 
"T 6 平整 冰 与 锥 体 结构 的 相互 作用 
e. 6 Interaction between level ice and conical structure 
Y sud h=0.15 0.5 m/ 
1-0.15 m v0.5 m/s 
a a H=0.65 m D=8.1 m 
GO z" 
N 
e 
CN 
> 
» 10 20 30 40 50 60 
ed 时 间 /s 
c 图 7 离散 元 模拟 锥 体 结 构 冰 载荷 时 程 
T Fig.7 Time series of ice load on conical structure 
(S) simulated with DEM 


为 验证 该 离散 元 计算 结果 的 合理 性 ,这 里 选取 
平整 冰 与 上 锥 风电 基础 结构 4 个 不 同 作用 部 位 的 冰 
载荷 峰值 的 均值 与 IEC 61400-3 规范 2 和 ISO 
19906 PRUE?! 进行 对 比 ,对 比 结果 如 图 8 所 示 。 从 
中 可 以 发 现 , 离 散 元 模拟 冰 载 荷 峰值 和 IEC 规范 及 
ISO 19906 标准 计算 冰 力 均 随 锥 径 的 增加 而 线性 增 
大 ,但 离散 元 结果 要 低 于 两 个 规范 值 , 其 中 ISO 
19906 标准 计算 的 冰 载 荷 最 大 。 这 是 由 于 ISO 标准 
和 IEC 规范 的 依据 主要 是 理论 模型 分 析 , 理论 模型 
一 般 不 考虑 动态 过 程 ,只 考虑 海 冰 断裂 瞬间 的 静态 
冰 力 ;也 不 考虑 圆锥 形状 ,只 考虑 二 维 的 斜坡 结构 。 
因此 理论 模型 对 冰 力 的 计算 一 般 较 大 。 此 外 ,考虑 
到 TEC 规范 及 ISO 19906 标准 相对 比较 保守 , 以 保 
障 结构 抗 冰 设计 的 安全 性 。 


ChinaXiv 合 作 期 刊 


9DEM ®IEC aISO 


20 
6.0 6.5 7.0 T3 8.0 8.5 
锥 径 /m 
8 ”离散 元 计算 冰 载 荷 与 ISO 标准 TEC 规范 对 比 
Fig.8 Comparison of ice loads of DEM results 
with ISO standard and IEC code 


2.2 ”潮汐 涨 落 影响 下 平整 冰 对 锥 体 结构 冰 载 荷 的 


影响 


受潮 汐 涨 落 影 响 , 海 冰 将 作用 于 上 下 不 同 锥 面 。 
为 分 析 平 整 冰 与 上 锥 、 下 锥 作用 下 破坏 模式 和 冰 载 荷 
的 差别 ,分 别 选取 水 位 高 度 太 为 0.65 m fil - 0.65 m, 
海 冰 作用 上 下 锥 径 均 为 D=8.1m。 图 9 给 出 了 平 
整 冰 对 风电 上 锥 和 下 锥 结构 相互 作用 的 离散 元 模拟 
过 程 , 计 算得 到 冰 载 荷 时 程 曲线 如 图 10 所 示 。 其 
中 ,上 锥 体 冰 载荷 峰值 均值 为 Fu =83.5 KN ,下 锥 体 
冰 载 集 峰 值 的 均值 为 Fu, =79. 1 KN。 不 难 发 现 ,不 
论 是 上 锥 体 还 是 下 锥 体 ,其 冰 载 荷 均 具有 很 强 的 随 
机 性 和 周期 脉冲 特点 。 冰 载荷 的 随机 性 是 由 海 冰 作 
用 于 锥 体 结构 上 时 发 生 的 随机 断裂 造成 的 , 其 周期 
性 与 冰 速 、 冰 厚 和 海 冰 强度 影响 下 的 断裂 长 度 密 不 
可 分 。 对 于 平整 冰 作 用 于 上 下 锥 时 冰 载 荷 的 差别 ， 
龙 雪 等 ”从 冰 排 断裂 时 的 受 力 特点 进行 了 详细 的 讨 
论 。 这 主要 是 由 于 海 冰 与 上 锥 体 作 用 时 会 发 生 明 显 
的 弯曲 破碎 , 破 雁 后 的 冰 块 在 后 续 冰 排 推动 下 不 断 
的 升 堆积 于 上 锥 面 ,后 受 重力 影响 落 入 水 中 被 清除 
的 周期 性 加 载 与 卸载 过 程 ;而 作用 于 下 锥 体 的 海 冰 
受 锥 体 结构 阻碍 而 不 断 下 潜 , 在 浮力 与 拖 忠 力 影 响 
下 , 碎 冰 块 对 锥 体 结构 的 作用 力 被 大 大 削弱 ,并 在 后 
续 海 冰 推 动 下 被 清除 ,完成 下 锥 冰 载 向 的 周期 性 加 
载 与 卸载 过 程 。 如 图 10 所 示 ,由 于 平整 冰 与 上 锥 持 
续 作 用 时 ,堆积 海 冰 未 能 及 时 清除 并 对 后 续 平 整 冰 
产生 运动 约束 , 相 较 于 下 锥 体 ,后 续 冰 力 不 能 快速 释 
放 , 从 而 导致 上 锥 冰 载 荷 高 于 下 锥 。 

图 11 给 出 了 风机 基础 结构 位 移 响 应 曲线 ,其 中 
上 锥 与 下 锥 位 移 响应 表现 出 相似 的 振动 模式 。 这 里 
截取 4 =25 ~45 s 时 的 风机 振动 状态 下 的 位 移 响应 ， 
上 锥 振动 的 位 移 平均 幅 值 4,, = 28.43 mm, 下 锥 振动 
的 位 移 平均 幅 值 4 =22. 92 mm。 此 外 ,不 难 发 现 


上 锥 仅 有 少 部 分 位 移 小 于 0, 而 下 锥 有 较 多 时 刻 出 
现 小 于 0 的 位 移 。 这 是 由 于 结构 运动 到 最 大 位 移 处 
冰 力 最 大 ,最 大 冰 力 作用 结束 后 冰 力 快速 释放 , 达到 
最 小 值 , 完 成 一 个 冰 力 周期 ,结构 迅速 回 弹 。 同 时 ， 
由 于 上 锥 最 低 冰 载荷 远 高 于 下 锥 ,使 得 结构 位 移 恢 
复 困难 ,难以 回 弹 至 运动 周期 的 初始 点 ,但 对 于 下 
锥 ,最 大 冰 载 荷 作 用 结束 后 冰 力 迅速 释放 接近 于 0, 
结构 相对 容易 回 弹 经 过 初始 点 。 但 无 论 如 何 , 海 冰 
作用 于 上 锥 时 的 振动 幅 值 要 大 于 下 锥 ,这 与 上 锥 结 
FK Bn] C T SC o 


) 平整 冰 与 上 锥 结构 相互 作用 (b) 平整 冰 与 下 锥 结构 相互 作用 


9 平整 冰 与 单 桩 锥 体 风 电 基 础 结构 的 离散 元 模拟 
ig. 9 DEM simulation of interaction between level ice and 
onical mono-pile structure of wind turbine foundation 
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h-0.15 m v=0.5 m/s 
H-0.65 m D-8.1 m 
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(a) 上 锥 体 冰 载荷 


chinaXiv 
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h=0.15m v=0.5 m/s 
H=-0.65 m D-8.1 m 


0 10 20 30 40 50 
时 间 /s 
(b) 下 锥 体 冰 载荷 
图 10 平整 冰 对 风电 基础 结构 冰 载 荷 的 离散 元 模拟 结果 


Fig. 10 Ice loads on wind turbine foundation structure 


simulated with DEM of level ice 
为 分 析 锥 径 对 风电 基础 结构 冰 载 荷 及 位 移 响应 
的 影响 ,选取 水 位 五 分 别 为 上 2.65 m、 +1.95 m, 
+1.30 m 和 + 上 0. 65 m 进行 海 冰 与 锥 体 结构 的 离散 元 
分 析 , 其 所 对 应 的 锥 径 D 为 6.15 m,6. 8 m,7.45 m 
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和 8.1m。 由 于 冰 和 载荷 时 程 的 随机 性 特点 ,本 文选 取 
结构 冰 载 丛 峰 值 的 均值 进行 对 比 ,如 图 12 所 示 。 结 
果 表明 ,相同 锥 径 下 上 锥 冰 和 载荷 明 显 大 于 下 锥 体 。 
在 较 大 锥 径 处 ,由 于 海 冰 作 用 位 置 处 的 锥 体 曲率 较 
小 , 易 引 发 海 冰 大 面积 拉 伸 弯曲 破碎 ,使 得 碎 冰 堆积 
体积 增加 ;而 在 较 小 锥 径 处 ,由 于 破碎 海 冰 较 少 , 海 
冰 作 用 位 置 处 的 锥 体 表面 曲率 较 大 ,摩擦 力 不 足 , 破 
碎 海 冰 易 被 清除 ,最 终 导 致 结构 冰 载 丛 随 锥 径 的 减 
小 而 减 小 1。 
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(b) 平整 冰 作 用 于 下 锥 体 
图 11 平整 冰 作 用 下 风电 基础 结构 的 位 移 响应 
Fig. 11 Displacement response of wind turbine foundation 
structure under action of level ice 
+ LHE 
85H 9 下 锥 


冰 载 荷 /kN 
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图 12 “不同 锥 径 下 平整 冰 对 风 构 结构 的 冰 载 荷 
Fig. 12 Ice loads of level ice on conical structure of 
wind turbine under various cone diameters 
图 13 给 出 了 平整 冰 作 用 下 锥 体 风电 基础 结构 
的 位 移 响应 随 锥 径 的 变化 。 可 以 发 现 ,无 论 海 冰 作 
用 于 上 锥 还 是 下 锥 , 风机 结构 的 位 移 响 均 随 锥 径 的 
增加 而 增 大 , 且 作 用 于 上 锥 时 的 振动 位 移 要 大 于 下 
锥 ,其 与 结构 冰 载 从 的 变化 规律 相似 。 
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图 13 锥 径 影响 下 风机 结构 的 冰 激 振动 位 移 
Fig. 13 The ice-induced displacements of structure vibration 


of wind turbine under influence of cone diameter 


3 和 碎 冰 区 锥 体 风 电 结构 冰 载荷 的 离散 
元 分 析 及 主要 影响 因素 分 析 


三 在 碎 冰 区 , 碎 冰 同 样 会 对 风电 的 振动 响应 与 作 
业 妥 全 产生 影响 。 为 研究 碎 冰 对 锥 体 结构 冰 载荷 的 
分 第 特性 ,对 不 同 密集 度 下 的 碎 冰 区 进行 离散 元 构 
造 3 计 算 分 析 碎 冰 作用 于 不 同 锥 体 部 位 处 的 结构 冰 
戎 符 与 结构 位 移 响应 的 变化 规律 。 


CD 
305 碎 冰 区 锥 体 风电 结构 冰 载 荷 的 离散 元 模拟 


——— a 

海 冰 作 用 于 上 锥 体 和 下 锥 体 情况 ,取水 位 高 度 
H'- 0. 65 m, 其 所 对 应 的 水 线 处 锥 径 D =8.1 m。 采 
用 移 展 多 面体 离散 元 方法 对 碎 冰 与 上 锥 和 下 锥 结 松 
HP e EE T EA n, 如 图 14 所 示 , 计 算 冰 
ed 
的 扰动 响应 如 图 16 所 示 。 


(a) 碎 冰 与 上 锥 体 相 互 作用 


(b) 碎 冰 与 下 锥 体 相 互 作 用 
图 14 碎 冰 与 单 桩 锥 体 风电 基础 结构 
相互 作用 的 离散 元 模拟 


Fig.14 Discrete element simulation of interactions between 


broken ice floes and conical mono-pile structure 


of wind turbine foundation 
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(b) 下 锥 体 冰 载荷 
图 15 碎 冰 对 风电 基础 结构 冰 载 荷 的 离散 元 模拟 
Fig.15 Ice loads on wind turbine structure simulated 


with DEM in broken ice field 
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(b) 雁 冰 区 下 锥 作用 位 移 时 程 
图 16 碎 冰 作用 下 风电 基础 结构 的 位 移 响应 
Fig.16 Displacement response of wind turbine foundation 

structure under action of broken ice floes 
碎 冰 区 结构 冰 载 倚 具 有 与 平整 冰 作 用 下 冰 和 载荷 
类 似 的 随机 性 特征 ,但 冰 载 荷 要 明显 低 于 平整 冰 ; 在 
碎 冰 区 ,上 锥 体 冰 载荷 峰值 的 均值 为 Fu =4.6 kN, 
而 下 锥 体 上 的 冰 载 荷 Fawn =5.6 KN ,明显 低 于 平整 
冰 作 用 下 的 已, =83.5 kN M Fann =79. 1 KN; REYK 
作用 于 上 锥 体 时 结构 振动 幅 值 的 平均 值 4,, = 


2.68 mm ,而 作用 于 下 锥 体 上 的 平均 值 Aion 22.82 mm, 
显著 低 于 平整 冰 作 用 下 的 Ap 228.43 mm 和 A, = 
22.92 mm。 此 外 ,还 可 以 发 现 ,在 碎 冰 区 作用 于 上 锥 
体 上 的 冰 载 荷 和 结构 振动 响应 均 低 于 下 锥 体 , 与 平 
整 冰 区 的 情况 正好 相反 。 这 主要 是 由 于 平整 冰 与 狼 
体 作 用 时 会 发 生 断 裂 ,上 锥 和 下 锥 上 破碎 冰 产 生 约 
束 条 件 的 不 同 导 致 了 上 锥 体 的 冰 载 荷 要 大 于 下 锥 
体 。 然 而 对 于 碎 冰 区 , 海 冰 不 会 发 生 断 裂 现象 , 碎 冰 
块 与 上 锥 作用 时 , 海 冰 在 后 续 不 连续 碎 冰 的 作用 下 
难以 向 上 人 攀 怜 而 在 海水 拖 遇 力作 用 下 更 容易 扩散 开 
来 ,导致 结构 冰 力 较 小 ; 碎 冰 块 与 下 锥 作用 时 碎 冰 下 
潜 , 相 较 于 上 锥 , 受 浮力 \ 拖 忠 力 的 影响 ,增加 了 碎 冰 
与 下 锥 的 连续 接触 作用 ,其 主要 表现 在 印 载 冰 力 较 
上 锥 释放 缓慢 ,进而 增加 了 下 锥 上 的 冰 力 。 因 此 , 碎 
水 的 上 能 冰 力 要 小 于 下 维 冰 力 。 


P 碎 冰 区 锥 径 及 海 冰 密 集 度 对 风电 
J lr] 


= HITER DX EEIE RT DUEB, Aul i 3) iK zB far Js 
EE ,分 别 选取 水 位 五 为 +2.60 m, +2.17 m, 
E 5m, £1.52 m, £1.30 m, £0.87 m fll £0.65 m 
时 进行 海 冰 与 锥 体 结构 的 离散 元 分 析 , 其 对 应 锥 径 
Dk 6. 15 m,6.5 m,6.8 m. 7.1 m,7. 45 m,7. 8 m 和 
8. 1m, 图 17 28:8. T ERYK yk ER far AE (ELT 35] (E DG E 

径 锡 变化 情况 。 图 18 给 出 了 碎 冰 区 锥 体 风 电 基 础 

结 移 位 移 响 应 随 锥 径 的 变化 情况 。 

6.5 


结构 冰 载 荷 


HE1 径 /m 
图 17 碎 冰 区 锥 径 对 结构 冰 载 荷 的 影响 


Fig.17 Influence of cone diameter on ice load 


in the broken ice field 

由 图 17 ~ 图 18 可 以 看 出 : 冰 载 荷 随 着 锥 径 的 增 
加 而 增 大 ,下 锥 冰 载 荷 大 于 上 锥 结构 。 这 是 因为 碎 
冰 与 上 锥 体 相互 作用 时 , 碎 冰 块 未 能 在 上 锥 面 发 生 
的 升 而 沿 锥 径 四 周 扩 散 开 来 ,如 图 14(a) 所 示 ; 248 
冰 与 下 锥 体 相互 作用 时 , 碎 冰 块 在 下 锥 面 的 阻碍 作 
用 下 发 生 下 潜 , 并 在 浮力 与 重力 的 共同 影响 下 使 其 
与 下 锥 面 持续 接触 作用 ,导致 水 力 增加 ,如 图 14(b) 
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所 示 。 因 此 , 碎 冰 对 下 锥 体 的 冰 载 荷 要 高 于 上 锥 。 
类 似 于 冰 载 荷 变化 规律 , 碎 冰 作 用 下 锥 结构 位 移 响 
应 要 大 于 上 锥 ,但 由 于 基础 结构 的 运动 刚度 较 大 , 导 
致 上 下 锥 运动 位 移 响应 的 差异 较 小 。 


3.0r 
LER: 
2.6| 9 F4 


ma 


位 移 /mm 


L L L L L 
6.0 6.5 7.0 7.5 8.0 8.5 
HEfS /m 


图 18 碎 冰 区 锥 径 对 结构 振动 位 移 的 影响 
Fig.18 Influence of cone diameter on vibration displacement 
of structure in the broken ice field 
为 研究 海 冰 密 集 度 对 上 锥 风电 结构 冰 载 荷 的 影 
啊 , 应 用 离散 元 方法 模拟 不 同 海 冰 密集 度 C 为 
6096 .65% 70% .75% 80% 85% 和 9096 EF ÉY YK ÆR 
fer IRRE 19 所 示 。 从 中 可 以 看 出 , 随 着 海 
冰 密 集 度 的 提高 ,结构 冰 载 荷 增加 显著 。 这 是 因为 
冰 块 在 较 低 密 集 度 下 与 结构 相互 作用 时 缺乏 有 效 约 
束 ,从 而 更 容易 从 锥 体 四 周 扩散 ,不 能 与 锥 体 结构 发 
生 连 续 作 用 ,从 而 降低 了 结构 冰 载 荷 。 随 着 海 冰 密 
集 度 的 增加 , 冰 块 间 的 相互 约束 不 断 增 大 ,从 而 也 增 
加 了 结构 冰 载 答 。 
7 


h=0.15m D=8.1 m 
6F v=0.5 m/s 4=2 m Ld 


50 60 70 80 90 100 
密集 度 /% 


图 19 海 冰 密集 度 对 冰 载 荷 的 影响 


Fig. 19 Influence of ice concentration on ice load 


4 平整 冰 和 碎 冰 区 单 桩 锥 体 风 电 基 础 
结构 冰 载 符 的 对 比分 析 


综合 以 上 离散 元 模拟 计算 结果 , 海 冰 对 风电 基 
结构 在 不 同 锥 径 下 的 冰 载 荷 和 结构 振动 响应 存在 
差异 。 图 20 给 出 了 离散 元 模拟 的 平整 冰 和 碎 冰 区 
上 下 锥 体 结构 冰 载 荷 情况 ,其 中 碎 冰 密集 度 C = 
80% ,平均 尺寸 4=2 m 。 不 难看 出 , 冰 载 位 均 随 锥 
径 的 增加 而 增加 。 平 整 冰 对 上 锥 结构 冰 载 荷 要 明显 
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大 于 对 下 锥 ,而 碎 冰 对 上 锥 结构 冰 载 荷 则 要 低 于 下 
锥 。 以 上 现象 产生 的 原因 在 上 文 进 行 了 讨论 说 明 。 
从 冰 载 人 答 的 幅 值 上 看 ,以 上 算 例 中 的 平整 冰 对 锥 体 
结构 的 冰 载 荷 是 碎 冰 区 的 10 倍 左右 。 当 然 , 浮 冰 块 
的 大 小 .密集 度 等 参数 对 冰 载 和 荷 的 影响 非常 显著 。 
因此 ,平整 冰 对 锥 体 结构 的 冰 载 荷 要 远 远 大 于 碎 冰 
区 ,其 主要 受 海 冰 强度 的 影响 。 为 此 ,在 对 冰 区 风电 
基础 的 结构 设计 和 安全 运行 , 则 需要 更 多 地 考虑 平 


整 冰 的 影响 。 
100 
90 
Ž so z 
3E 70 E 
bu E 
关 60 z 
Xx 办 
> 30 3 
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C20 不 同 锥 径 下 平整 冰 和 碎 冰 对 锥 体 结构 的 冰 载 茶 


-Fig.20 Ice loads of conical structure of level ice and 


= broken ice floes under various cone diameters 
EO 入 分 析 在 不 同 难 径 条 件 下 海 冰 对 锥 体 风 电 基础 
结 峰 振动 响应 的 影响 ,对 两 种 冰 况 作用 下 的 结构 位 
移 测 应 进行 对 比分 析 , 如 图 21 所 示 。 可 以 发 现 ,在 
礁 同 冰 况 条 件 下 ,平整 冰 区 结构 位 移 幅 值 响应 明显 
高 至 碎 冰 区 ;而 在 相同 冰 况 条 件 下 ,结构 位 移 随 锥 径 


35 7 
Pr 30r a 下 锥 (平整 冰 ) 6 
E al e 上 锥 ( 碎 冰 ) E 
O 二 | o FEK E 
E» ig 
SE 3 
X = 
88 10 2 
E u 
5 1 
0 0 
3/5 6.0 6.5 7.0 T5 8.0 8.5 
锥 径 /m 
了 21 ”不同 锥 径 下 平整 冰 和 碎 冰 作用 下 锥 体 结构 的 振动 位 移 
Fig.21 Vibration displacements of conical structure of level 


ice and broken ice under various cone diameters 
Zt VN 
5 £i it 
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散 元 方法 计算 分 析 了 单 桩 锥 体 风 电 基 础 结构 分 别 与 
平整 水 和 碎 冰 的 相互 作用 过 程 ,并 通过 与 IEC 规范 
和 ISO 标准 对 计算 的 冰 载 荷 进 行 了 对 比 验证 。 采 用 
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离散 元 方法 分 析 了 平整 冰 与 碎 冰 对 锥 体 不 同位 置 处 
的 冰 和 载荷 和 冰 激 结 构 响 应 ,得 到 以 下 结论 。 

1) 冰 载荷 随 作 用 锥 径 的 增加 而 增 大 ,平整 冰 对 
锥 体 结构 的 冰 载 荷 要 明显 大 于 碎 冰 。 
2) 平 整 冰 对 上 锥 结构 的 冰 和 载荷 要 明显 高 于 下 
锥 ,而 雁 冰 对 上 锥 作用 冰 载 荷 要 低 于 下 锥 。 

3) 扩 展 多 面体 离散 元 方法 可 有 效 模 拟 海 冰 对 单 
桩 锥 体 风电 结构 的 相互 作用 过 程 ,准确 计算 分 析 冰 
载荷 的 变化 规律 ,可 为 冰 区 锥 体 风 电 结 构 的 抗 冰 设 
计 提 供 有 益 参 考 。 
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